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Abstract 

Exclusive hadron production processes in photon-photon collisions provide impor- 
tant tests of QCD at the amplitude level, particularly as measures of hadron distri- 
bution amplitudes and skewed parton distributions. The determination of the shape 
and normalization of the distribution amplitudes has become particularly important 
in view of their importance in the analysis of exclusive semi-leptonic and two-body 
hadronic 5-decays. Interesting two-photon physics, including doubly-tagged 7*7* 
reactions, will be accessible at low energy, high luminosity e^e~ colliders, including 
measurements of channels important in the light-by- light contribution to the muon 
g~2 and the study of the transition between threshold production controlled by low- 
energy effective chiral theories and the domain where leading-twist perturbative QCD 
becomes applicable. The threshold regime of hadron production in photon-photon 
and e^e~ annihilation, where hadrons are formed at small relative velocity, is par- 
ticularly interesting as a test of low energy theorems, soliton models, and new types 
of resonance production. Such studies will be particularly valuable in double-tagged 
reactions where polarization correlations, as well as the photon virtuality dependence, 
can be studied. 
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1 Introduction 

Two-photon annihilation 7*(q'i)7*(5'2) hadrons for real and virtual photons provide 
some of the most detailed and incisive tests of QCD. Among the processes of special 
interest are: 

1. the total two-photon annihilation hadronic cross section cr(s,g^,g|), which is 
related to the light-by-light hadronic contribution to the muon anomalous mo- 
ment; 

2. the formation of C = -|- hadronic resonances, which can reveal exotic states 
such as qqg hybrids and discriminate gluonium formation 0; 

3. single-hadron processes such as 7*7* vr°, which test the transition from the 
anomaly-dominated pion decay constant to the short-distance structure of cur- 
rents dictated by the operator-product expansion and perturbative QCD fac- 
torization theorems; 

4. hadron pair production processes such as 7*7 n^n^ , K^K~ ,pp, which at 
fixed invariant pair mass measures the s — > t crossing of the virtual Compton 
amplitude [§, When one photon is highly virtual, these exclusive hadron 
production channels are dual to the photon structure function F2{x,Q'^) in 
the endpoint x —>■ 1 region at fixed invariant pair mass. The leading twist- 
amplitude for 7*7 — i> '7r"'"7r~ is sensitive to the — 1/(1— x) moment of 
the qq distribution amplitude of the two-pion system p, ||, the 
timelike extension of skewed parton distributions. In addition one can measure 
the pion charge asymmetry in e~^e~ 7r^7r~e^e~ arising from the interference 
of the 77 TT+TT^ Compton amplitude with the timelike pion form factor |^. 
At the unphysical point s = qf = q2 = 0, the amplitude is fixed by the low 
energy theorem to the hadron charge squared. As reviewed by Karliner in these 
proceedings 0] , the ratio of the measured 77 AA and 77 — > pp cross sections 
is anomalous at threshold, a fact which may be associated with the soliton 
structure of baryons in QCD 0; 

5. At large momentum transfer, the angular distribution of hadron pairs produced 
by photon-photon annihilation are among the best determinants of the shape of 
the meson and baryon distribution amplitudes (j)M{x,Q), and ^^(xj,^) which 
control almost all exclusive processes involving a hard scale Q. The determi- 
nation of the shape and normalization of the distribution amplitudes, which 
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are gauge-invariant and process-independent measures of the valence wavefunc- 
tions of the hadrons, has become particularly important in view of their impor- 
tance in the analysis of exclusive semi-leptonic and two-body hadronic i?-decays 
14, 15 1 . There has also been considerable progress both in cal- 



10, 11, 13, \U 



culating hadron wavefunctions from first principles in QCD and in measuring 
them using diffractive di-jet dissociation. 



Much of this important two-photon physics is accessible at low energy, high lumi- 
nosity e^e~ colliders such as the proposed PEP-N project, particularly for measure- 
ments of channels important in the light-by-light contribution to the muon g-2 and 
the exploration of the transition between threshold amplitudes which are controlled 
by low-energy effective theories such as the chiral Hamiltonian through the transition 
to the domain where leading-twist perturbative QCD becomes applicable. There have 
been almost no measurements of double-tagged events needed to unravel the sepa- 
rate qf and g| dependence of photon-photon annihilation. Hadron pair production 
from two-photon annihilation plays a crucial role in unraveling the perturbative and 
non-perturbative structure of QCD, first by testing the validity and empirical appli- 
cability of leading-twist factorization theorems, second by verifying the structure of 
the underlying perturbative QCD subprocesses, and third, through measurements of 
angular distributions and ratios which are sensitive to the shape of the distribution 
amplitudes. In effect, photon-photon collisions provide a microscope for testing fun- 
damental scaling laws of PQCD and for measuring distribution amplitudes. It would 
also be interesting to measure the novel relativistic atomic coalescence processes, sin- 
gle and double muonium formation: e'^e~ — > [/i"'"e~]/i+e~ and e+e~ — > [//"""e"] [//"""e"] . 



2 The Photon-to-Pion Transition Form Factor and 
the Pion Distribution Amphtude 

The simplest and perhaps most elegant illustration of an exclusive reaction in QCD is 
the evaluation of the photon-to-pion transition form factor F^^t^{Q'^) which is measur- 
able in single-tagged two-photon ee — > eevr^ reactions. The form factor is defined via 
the invariant amplitude F'^ = —ie^ F.,^^{Q'^)e^"P'' ple pq^j . As in inclusive reactions, one 
must specify a factorization scheme which divides the integration regions of the loop 
integrals into hard and soft momenta, compared to the resolution scale Q. At leading 
twist, the transition form factor then factorizes as a convolution of the 7*7 qq 
amplitude (where the quarks are coUinear with the final state pion) with the valence 
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light-cone wavefunction of the pion [^: 



FMQ') = ^ I' dxMx, Q)T%,{x, g2). (1) 

The hard scattering amphtude for 77* qq is T^^j^{x, Q^) = [(1 — x)Q'^]~'^ (1 + 0{as))- 
For the asymptotic distribution amphtude cj)'^^'^^^ [x] = ^/SfjrX^l — x) one predicts 
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Q^F^^iQ^) = 2/J 1 



where Q* = e ^^^Q is the estimated BLM scale for the pion form factor in the V 
scheme. 

The PQCD predictions have been tested in measurements of 67 evr" by the 
CLEO collaboration ^ (see Figure |l| (b)). The flat scaling of the Q^F^^{Q'^) data 



from Q"^ = 2 to = 8 GeV^ provides an important confirmation of the applicability 
of leading twist QCD to this process. The magnitude of is remarkably 

consistent with the predicted form, assuming the asymptotic distribution amplitude 
and including the LO QCD radiative correction with av{e~^^'^Q) /tt ~ 0.12. One could 
allow for some broadening of the distribution amplitude with a corresponding increase 
in the value of ay at small scales. Radyushkin [|l8l , Ong |jl9[ and Kroll have also 
noted that the scaling and normalization of the photon to pion transition form factor 
tends to favor the asymptotic form for the pion distribution amplitude and rules out 



broader distributions such as the two- humped form suggested by QCD sum rules pl j. 
More comprehensive analyses, which include consideration of next-to-leading order 
corrections and some higher-twist contributions dictated by vector meson spectra 
and QCD sum rules have been given by A. Khodjamirian ||22|, A. Schmedding and 
O. Yakovlev [H, and by A. P. Bakulev et al. 



When both photons are virtual, the denominator of Th for the 77* — > tt*' reaction 
becomes (1 — x)Q\ + XQ2 @, [T^) and the amplitude becomes nearly insensitive to the 
shape of the distribution amplitude once it is normalized to the pion decay constant. 
Thus the ratio of singly virtual to doubly virtual pion production is particularly 
sensitive to the shape of (j).jr{x,Q'^) since higher order corrections and normalization 
errors tend to cancel in the ratio. 
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Figure 1: (a) Transverse lattice results for the pion distribution amplitude at 
g2 ^ lOGeVl The solid curve is the theoretical prediction from the combined 
DLCQ/transverse lattice method the chain line is the experimental result ob- 
tained from dijet diffractive dissociation pBl E^. Both are normalized to the same 



area for comparison, (b) Scaling of the transition photon to pion transition form fac- 
tor Q'^F^^o{Q^). The dotted and solid theoretical curves are the perturbative QCD 
prediction at leading and next-to-leading order, respectively, assuming the asymptotic 
pion distribution The data are from the CLEO collaboration []T7 |. 
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3 Non-Perturbative Calculations of the Pion Dis- 
tribution Amplitude 

The distribution amplitude (f){x, Q) can be computed from the integral over transverse 
momenta of the renormalized hadron valence wavefunction in the light-cone gauge at 
fixed light-cone time 0]: 

<P{x,Q) = J dPk^e {q' - ^^^^ i^^'^KxM). (2) 

where a global cutoff in invariant mass is identified with the resolution Q. The distri- 
bution amplitude Q) is boost and gauge invariant and evolves in hiQ through an 
evolution equation [Q. Since it is formed from the same product of operators as the 
non-singlet structure function, the anomalous dimensions controlling Q) depen- 
dence in the ultraviolet log Q scale are the same as those which appear in the DGLAP 
evolution of structure functions [^. The decay n ^ fiu normalizes the wave func- 
tion at the origin: dx(j){x , Q) = /^/(2^/3). One can also compute the distribution 
amplitude from the gauge invariant Bethe-Salpeter wavefunction at equal light-cone 
time. This also allows contact with both QCD sum rules and lattice gauge the- 
ory; for example, moments of the pion distribution amplitudes have been computed 
in lattice gauge theory Conformal symmetry can be used as a template 



to organize the renormalization scales and evolution of QCD predictions [pSj |33| . For 
example, Braun and collaborators have shown how one can use conformal symmetry 
to classify the eigensolutions of the baryon distribution amplitude ||3^ . 



Dalley and Burkardt and Seal have calculated the pion distribution am- 



plitude from QCD using a combination of the discretized light-cone quantization 
method for the x~ and x'^ light-cone coordinates with the transverse lattice method 
P7|, RBI in the transverse directions, A finite lattice spacing a can be used by choosing 



the parameters of the effective theory in a region of renormalization group stability to 
respect the required gauge, Poincare, chiral, and continuum symmetries. The overall 
normalization gives /^r = 101 MeV compared with the experimental value of 93 MeV. 
Figure |I| (a) compares the resulting DLCQ/transverse lattice pion wavefunction with 
the best fit to the diffractive di-jet data (see the next section) after corrections for 
hadronization and experimental acceptance The theoretical curve is somewhat 



broader than the experimental result. However, there are experimental uncertainties 
from hadronization and theoretical errors introduced from finite DLCQ resolution, 
using a nearly massless pion, ambiguities in setting the factorization scale Q^, as well 
as errors in the evolution of the distribution amplitude from 1 to 10 GeV^. Instan- 
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ton models also predict a pion distribution amplitude close to the asymptotic form 
p9| . In contrast, recent lattice results from Del Debbio et al. predict a much 
narrower shape for the pion distribution amplitude than the distribution predicted 
by the transverse lattice. A new result for the proton distribution amplitude treating 
nucleons as chiral solitons has recently been derived by Diakonov and Petrov ^ . 
Dyson- Schwinger models of hadronic Bethe-Salpeter wavefunctions can also be 
used to predict light-cone wavefunctions and hadron distribution amplitudes by in- 
tegrating over the relative k~ momentum. There is also the possibility of deriving 
Bethe-Salpeter wavefunctions within light-cone gauge quantized QCD [^] in order 
to properly match to the light-cone gauge Fock state decomposition. 



4 Measurements of the Pion Distribution Ampli- 
tude by Di-jet Diffract ive Dissociation 

The shape of hadron distribution amplitudes can be measured in the diffractive dis- 
sociation of high energy hadrons into jets on a nucleus. For example, consider the 



reaction P3| , ^ ^ tiA Jeti + Jet2 + A' at high energy where the nucleus A' is 
left intact in its ground state. The transverse momenta of the jets balance so that 
+ A;_L2 = if± < R~^A- The light-cone longitudinal momentum fractions also need 
to add to Xi + X2 ~ 1 so that Api <R~A^. The process can then occur coherently in 
the nucleus. Because of color transparency and the long coherence length, a valence 
qq fluctuation of the pion with small impact separation will penetrate the nucleus 
with minimal interactions, diffracting into jet pairs |^3|]. The Xi = x, X2 = 1 — x 



dependence of the di-jet distributions will thus reflect the shape of the pion valence 
light-cone wavefunction in x; similarly, the k^i — k±2 relative transverse momenta of 
the jets gives key information on the second derivative of the underlying shape of the 
valence pion wavefunction ^ ^ . The diffractive nuclear amplitude extrapolated 



to t = should be linear in nuclear number A if color transparency is correct. The 
integrated diffractive rate should then scale as A^/R^ ~ A^^^. 

The E791 collaboration at Fermilab has recently measured the diffractive di-jet 



dissociation of 500 GeV incident pions on nuclear targets p6[. The results are consis- 
tent with color transparency, and the momentum partition of the jets conforms closely 
with the shape of the asymptotic distribution amplitude, 0^'^y™P*(x) = -\/3/^x(l — x), 
corresponding to the leading anomalous dimension solution to the perturbative 
QCD evolution equation. 

The interpretation of the diffractive dijet processes as measures of the hadron 
distribution amplitudes has recently been questioned by Braun et al. |^ and by 
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Chernyak who have calculated the hard scattering amplitude for such processes 
at next-to-leading order. However, these analyses neglect the integration over the 
transverse momentum of the valence quarks and thus miss the logarithmic ordering 
which is required for factorization of the distribution amplitude and color filtering in 
nuclear targets. 



5 Exclusive Two-Photon Annihilation into Hadron 
Pairs 




Figure 2: Comparison of the sum of 77 tt+tt^ and 77 K^K~ meson pair 
production cross sections with the scaling and angular distribution of the perturbative 
QCD prediction [Q]. The data are from the CLEO collaboration [pO|] . 

Two-photon reactions, 77 — HH at large s = {ki + A;2)^ and fixed Ocm-, pro- 
vide a particularly important laboratory for testing QCD since these cross-channel 
"Compton" processes are the simplest calculable large-angle exclusive hadronic scat- 
tering reactions. The helicity structure, and often even the absolute normalization 
can be rigorously computed for each two-photon channel [Q]. In the case of meson 
pairs, dimensional counting predicts that for large s, s'^da / dt{'y'y — >■ MM scales at 
fixed t/s OT 6cm up to factors of In s/ . The angular dependence of the 77 HH 
amplitudes can be used to determine the shape of the process-independent distribu- 
tion amplitudes, (j)H{x,Q). An important feature of the 77 — > MM amplitude for 
meson pairs is that the contributions of Landshoff pitch singularities are power-law 
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suppressed at the Born level — even before taking into account Sudakov form factor 
suppression. There are also no anomalous contributions from the a; ^ 1 endpoint 
integration region. Thus, as in the calculation of the meson form factors, each fixed- 
angle helicity amplitude can be written to leading order in l/Q in the factorized form 
[g2 = p'^ = tu/s; Qx = min(a;Q, (/ - x)Q)]: 

4>M{y, Qy)TH{x, y, S, 9cm4>M{x, Q^), (3) 

where Th is the hard-scattering amplitude 77 {,<i<i){,'i<i) for the production of the 
valence quarks collinear with each meson, and 4>m{x, Q) is the amplitude for finding 
the valence q and q with light-cone fractions of the meson's momentum, integrated 
over transverse momenta k± < Q. The contribution of non- valence Fock states are 



power-law suppressed. Furthermore, the helicity-selection rules [^] of perturbative 
QCD predict that vector mesons are produced with opposite helicities to leading order 
in 1/Q and all orders in a^. The dependence in x and y of several terms in Tx^y is 
quite similar to that appearing in the meson's electromagnetic form factor. Thus 
much of the dependence on (j)M{x,Q) can be eliminated by expressing it in terms 
of the meson form factor. In fact, the ratio of the 77 vr+vr" and e^e~ fi^fi~ 
amplitudes at large s and fixed 9cm is nearly insensitive to the running coupling and 
the shape of the pion distribution amplitude: 



s 



f (77 1-C0S2^, 



(4) 



cm 



The comparison of the PQCD prediction for the sum of tt+tt^ plus K^K~ channels 
with recent CLEO data 0] is shown in Figure |^. The CLEO data for charged 
pion and kaon pairs show a clear transition to the scaling and angular distribution 
predicted by PQCD for W = \J\s^^ > 2 GeV. It is clearly important to measure 
the magnitude and angular dependence of the two-photon production of neutral pions 
and p'^p~ in view of the strong sensitivity of these channels to the shape of meson 
distribution amplitudes (see Figures |]and H). QCD also predicts that the production 
cross section for charged p-pairs (with any helicity) is much larger than for that of 
neutral p pairs, particularly at large angles. Similar predictions are possible for 
other helicity-zero mesons. 

The analysis of exclusive B decays has much in common with the analysis of 
exclusive two-photon reactions For example, consider the three representative 
contributions to the decay of a. B meson to meson pairs illustrated in Figure |^. 
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Figure 3: Predictions for the angular distribution of the 77 tt+tt and 77 vr'^vr'' 
pair production cross sections for three different pion distribution amphtudes 

In Figure H(a) the weak interaction effective operator O produces a gg in a color 
octet state. A gluon with virtuality = 0{M'^) is needed to equilibrate the large 
momentum fraction carried by the b quark in the B wavefunction. The amplitude 
then factors into a hard QCD/electroweak subprocess amplitude for quarks which are 
collinear with their respective hadrons: T//([6(a;)n(l —x)] [q{y)u{l — y)]i[q{z)q{l — 
z)]2) convoluted with the distribution amplitudes Q) [^] of the incident and final 
hadrons: 

MoctetiB M1M2) = C dz dy dx 

Jo Jo Jo 

Q)Th{x, y, z)(t)MAy^ Q)<pM2{.z, Q). 

Here x = /pj^ = {k^ + / {p^h^P^i) light-cone momentum fractions carried 

by the valence quarks. 

There are a several features of QCD which are required to ensure the consistency of 
the PQCD approach: (a) the effective QCD coupling as{Q'^) needs to be under control 
at the relevant scales of B decay; (b) the distribution amplitudes of the hadrons need 
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z2 = C0S2 (9) 



Figure 4: Predictions for the angular distribution of the 77 p+p and 77 p^p^ 
pair production cross sections for three different p distribution amphtudes as in Figure 



to satisfy convergence properties at the endpoints; and (c) one requires the coherent 
cancelation of the couplings of soft gluons to color-singlet states. This property, 
color transparency |52[, is a fundamental coherence property of gauge theory and 
leads to diminished final-state interactions and corrections to the PQCD factorizable 
contributions. The problem of setting the renormalization scale of the coupling for 
exclusive amplitudes is discussed in |TB| . 



Baryon pair production in two-photon annihilation is also an important testing 
ground for QCD. The calculation of Th for Compton scattering requires the evaluation 
of 368 hehcity-conserving tree diagrams which contribute to 'y{qqq) — ^ l\qqq)' at the 
Born level and a careful integration over singular intermediate energy denominators 
p3| , 0, |55|. Brooks and Dixon |56| have recently completed a recalculation of the 



proton Compton process at leading order in PQCD, extending and correcting earlier 
work. It is useful to consider the ratio da/dt^jj pp)/da/dt{e~^e~ pp) since 
the power-law fall-off, the normalization of the valence wavef unctions, and much of 
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Figure 5: Three representative contributions to exclusive B decays to meson pairs in 
PQCD. The operators O represent the QCD-improved effective weak interaction. 



the uncertainty from the scale of the QCD coupling cancel. The scaling and angular 
dependence of this ratio is sensitive to the shape of the proton distribution amplitudes. 
The perturbative QCD predictions for the phase of the Compton amplitude phase can 
be tested in virtual Compton scattering by interference with Bethe-Heitler processes 

It is also interesting to measure baryon and isobar pair production in two photon 
reactions near threshold. Ratios such as (7(77 )/o'{'~f'-f A+A^) can be 

as large as 16 : 1 in the quark model since the three-quark wavefunction of the A 
is expected to be symmetric [||. Such large ratios would not be expected in soliton 
models M in which intermediate multi-pion channels play a major role. 



Recently Pobylitsa et al. have shown how the predictions of perturbative 



QCD can be extended to processes such as 77 — >• ppvr where the pion is produced at 
low velocities relative to that of the p or p by utilizing soft pion theorems in analogy 
to soft photon theorems in QED. The distribution amplitude of the pn composite is 
obtained from the proton distribution amplitude from a chiral rotation. A test of 
this procedure in inelastic electron scattering at large momentum transfer ep pix 
and small invariant p'n mass has been remarkably successful. Many tests of the soft 
meson procedure are possible in multiparticle e^e~ and 77 final states. 



6 Conclusions 

The leading-twist QCD predictions for exclusive two-photon processes such as the 
photon-to-pion transition form factor and 77 — hadron pairs are based on rigorous 
factorization theorems. The recent data from the CLEO collaboration on F^t^{Q'^) 
and the sum of 77 —>■ vr+vr" and 77 K^K~ channels are in excellent agreement 
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with the QCD predictions. It is particularly compelling to see a transition in angular 
dependence between the low energy chiral and PQCD regimes. The success of leading- 
twist perturbative QCD scaling for exclusive processes at presently experimentally 
accessible momentum transfer can be understood if the effective coupling av{Q*) is 
approximately constant at the relatively small scales Q* relevant to the hard scattering 
amplitudes [0. The evolution of the quark distribution amplitudes in the low-Q* 
domain also needs to be minimal. Sudakov suppression of the endpoint contributions 
is also strengthened if the coupling is frozen because of the exponentiation of a double 
logarithmic series. 

One of the formidable challenges in QCD is the calculation of non-perturbative 
wavefunctions of hadrons from first principles. The recent calculation of the pion 
distribution amplitude by Dalley |]25| and by Burkardt and Seal using light-cone 
and transverse lattice methods is particularly encouraging. The predicted form of 
07r(a:, Q) is somewhat broader than but not inconsistent with the asymptotic form 
favored by the measured normalization of Q'^F^^o{Q'^) and the pion wavefunction 
inferred from diffractive di-jet production. 

Clearly much more experimental input on hadron wavefunctions is needed, partic- 
ularly from measurements of two-photon exclusive reactions into meson and baryon 
pairs at the high luminosity B factories. For example, as shown in Figure 0, the ratio 



TT TT 



is particularly sensitive to the shape of pion distribution amplitude. At fixed pair 
mass, and high photon virtuality, one can study the distribution amplitude of multi- 
hadron states 0. Two-photon annihilation will provide much information on fun- 
damental QCD processes such as deeply virtual Compton scattering and large angle 
Compton scattering in the crossed channel. I have also emphasized the interrelation 
between the wavefunctions measured in two-photon collisions and the wavefunctions 
needed to study exclusive B and D decays. 

Much of the most interesting two-photon annihilation physics is accessible at low 
energy, high luminosity e^e~ colliders, including measurements of channels important 
in the light-by-light contribution to the muon g-2 and the study of the transition 
between threshold production controlled by low-energy effective chiral theories and 
the domain where leading-twist perturbative QCD becomes applicable. 

The threshold regime of hadron production in photon-photon and e~^e~ annihila- 
tion, where hadrons are formed at small relative velocity, is particularly interesting 
as a test of low energy theorems, soliton models, and new types of resonance pro- 
duction. Such studies will be particularly valuable in double-tagged reactions where 
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polarization correlations, as well as the photon virtuality dependence, can be studied. 
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